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Abstract. Contrary to what would be predicted on the basis of Cramer's model concerning 
the distribution of prime numbers, we develop evidence that the distribution of xp{x + H) — 
xp{x), for < a; < A^, is approximately normal with mean ~ H and variance ~ H log N/H, 
when < H < N^-\ 



0. Introduction 

Cramer [4] modeled the distribution of prime numbers by independent random variables 
Xn (for n > 3) that take the value 1 (n is "prime") with probability 1/logn and take 
the value {n is "composite") with probability 1 — 1/logn. If pn denotes the n^^ prime 
number this model predicts that 

lim ^card{n : 1 < n < N,pn+i - Pn > clogpn} = e"'' 

for all fixed positive real numbers c. Gallagher [6] showed that the above follows from 
Hardy & Littlewood's [10, p. 61] quantitative version of the prime /c-tuple conjecture: If 
D = {di, d2, dk} is a, set of k distinct integers, then 

k 

(1) Yl n^(^+^^) = i&m+o{i))x 

n<.x 1=1 

as a; — s> cxD where is the singular series 

- E (nfl) E KE^) 

qi,---,qk ai,...,afe i^i "i* 

l<i3'i<cxD l<ai<qi 

{ai,qi) = l 

E O'i/qiEi' 
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where e{9) = e^^^^. Hardy &: Littlewood showed that the right hand side above may be 
written more transparently as 



where i^p{D) denotes the number of distinct residue classes modulo p found among the 
members of D. (See the remarks following the proof of Lemma 3 in §2.) Since iyp{D) — k 
for all sufficiently large p, the product (3) is absolutely convergent. Hence its value is if 
and only if there is a prime p for which I'pi'D) = p. Gallagher [6] showed that from (1) it 
follows that 



when A X 1. Here m/c(A) = E(y^) is the /c*^ moment of a Poisson random variable Y with 
parameter A, and a x 6 means that a/b lies between two positive absolute constants. Thus 
the distribution of tt{x + h) — 7r{h) is approximately Poisson when h x logX, as predicted 
by the Cramer model. 

In this paper we investigate the distribution of primes in longer intervals. Let H = H{N) 
be a function of such that H — o{N) and H/ log N ^ oo as N ^ cxo. The Cramer model 
predicts that the distribution of 'ip{n + H) — ip{n) (for n < N) is approximately normal 
with mean ~ H and variance ~ HlogN. Assuming a strong form of the Hardy-Littlewood 
conjecture (1) we will show that this prediction holds in the range where H/logN oo 
and log if/ log — > as oo. In the range A^*^ < H < N^~^ we provide evidence 

showing that the Cramer model is incorrect, and conjecture instead that the distribution 
of 7/^(77- + H) —^tpin) is approximately normal with mean ~ H and variance ~ H \og{N/ H). 

When h x logx, the moments of ^ip^x + h) — il^{x) and of '^{x + h) — "(p^x) — h are of 
the same order of magnitude. However, for larger h one would expect that the moments 
of il){x + h) — i/j{x) to be far larger than those of i/j{x + h) — i/j{x) — h. We obtain our 
conclusions on the distribution of •0(a; + /i) — ^p{x) by analyzing these latter, more delicate 
moments. To facilitate this study, we set Ao(n) = A(n) — 1, with the result that 



Thus the main term is eliminated at the outset, which simplifies our calculations consider- 
ably. We recast (1) in an equivalent form that pertains to Ao(n): If di, . . . , are distinct 
integers, then 






{iPix + A log a;) - i;{x))'' dx mk{X)X{\ogX)'' 





x<n<x+h 



k 



(4) 




as a; — > 00 where &o{D) is related to 6(D) by the identities 




3CD 
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Here it is to be understood that 60 (0) = ©(0) = 1- One of the main steps in Gallagher's 
argument is to show that 

(7) ®(^) - 

di,... ,dk 
l<di<h 
di distinct 

as /i — > 00. There are (^) subsets D under consideration, but each one occurs k\ times in 
the above sum. Thus the above asserts that the mean value of 6(D) tends to 1 as /i — > 00. 
Correspondingly, we need to estimate the quantities 

(8) Rk{h)= 



di ,dk 
l<di<h 
dj distinct 



From (2) and (5) we see that 



qi,---,qk ai,...,ak i=l ^* 

l<gi<c>o l<ai<qi 

{ai,qi) = l 
E ailqi^l. 

The task of estimating averages of this expression is quite challenging, but our burden is 
substantially lightened by work of Montgomery & Vaughan [16] concerning a strikingly 
similar quantity. Let 



(10) mu{si;h) = Y\ E 1 - 



n=l \ ra=\ 

(rn+n,q) = \ 



be the k^^ centered moment of the number of reduced residues (mod q) in an interval. 
Lemma 2 of Montgomery & Vaughan asserts that 

(11) mfe(9;/.) = Q(^)Vfe(?;/.) 



where 



\<di<h l<9ik l<ai<qi 

(ai,Qi) = l 
E o,ilqi&L 

When k — 1, the conditions in the innermost sum cannot be fulfilled, and thus Vi{q; h) = 0. 
When k = 2, the conditions in the innermost sum require that = = oi + 02- Thus 

(13) ''^te'') = EsC t l«(''/<')l' 

d\q ' «=1 
d>l (a,d)=l 
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where 



(14) 



h 

E(a) = e{ma). 

m=l 



Montgomery & Vaughan showed that 



for each positive integer k. Unfortunately, this is not quite sharp enough for our present 
purposes, so our first job is to refine the above. 

Theorem 1. In the above notation, 



for every positive integer k, where //^ = 1 • 3 • • • (A; — 1) if k is even, and /ik = ^ if k is 
odd. 

Here the main term is the k^^ moment of a normal random variable with expectation 
and variance V2{q; h). We remark that the work of Granville & Soundararajan [7] (see 
§6a) places restrictions on the uniformity (in k) with which (15) can possibly hold. With 
Theorem 1 in hand, we are able to estimate the Rk{h). 

Theorem 2. Let h be an integer, h > 1, and suppose that Rk{h) is defined as in (8). 



for any nonnegative integer k, where A = 2 — Cq — log 2ti and Co denotes Euler's constant. 

For the smallest values of /c, one can be more precise, since it is clear that Ro{h) = 1, 
and that Ri{h) = 0. Also, from (5) and (48) we know that 



(15) 




Then 



Rk{h) = fXk{-hlogh + Ahfl'' + Ofc(/i^/2-V(7fe)+=) 



(16) 



R^ih) = -hlogh + Ah + 0{h^/'^+^). 



From (6) it follows that the left hand side of (7) is 




Hence we obtain Gallagher's estimate (7) in the more precise form 



(17) 




(2) (1 - Co - log27r)/i'=-i + 0(/i'=-W=). 



l<di<h 
dj distinct 
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We put 

(18) Mk{N; H) = Y^iHn + H) - V'(n) - i?)^, 



N 



and note that this is 



N H ^ N K 



(19) =E(E^o(n + /i)) = J] J] n^o(n + /i.)- 

n=l hi,...,hK n=l i=l 

l<hi<H 

Here the /i^ are not necessarily distinct, but once the distinct vahics have been identified, 
and their multiphcities accounted for, we can appeal to (4). The quantities Rk{h) arise in 
the main term. Thus from Theorem 2 we can derive an asymptotic estimate for the above, 
provided that the error term in (4) is sufficiently small and H is not too large. 

Theorem 3. Let Ek{x;'D) be defined by the relation 

k 

J2 l[Hn + di) = e{T>)x + Ek{x;'D), 

n<x i=l 

and suppose that 

(20) EkixiT)) <^ N^/^+^ 

uniformly for l<k<K,0<x<N, and distinct di satisfying 1 < di < H. Then 

r-TV 



{logx/ H + B)^/^dx 
(21) + O (iV(log iV)^/2if^/2 (^JL-^ '^'^^""^ + iy^Ari/2+.^ 



uniformly for log AT < H < N^/^ , where B = 1 — Cq — log27r and Cq denotes Euler's 
constant. 

In the case k = 1, the set D is a singleton, 6(23) = 1, and the hypothesis that 
Ei{N: {1}) ^ A^V2+e is equivalent to the Riemann Hypothesis (RH). 

In place of (20) if we assume only that Ei.{x; D) <C £^ for some E > Ari/2+£ then (21) 
holds with the modified error term 



We note that in deriving Theorem 3 from Theorem 2, we start with hi that are not 
necessarily distinct, and must reduce to distinct di, which is the reverse of the problem 
encountered in deriving Theorem 2 from Theorem 1, where we start with distinct di, and 
want to appeal to an estimate involving not necessarily distinct m^. 

Since the fiK are the moments of a normal random variable with mean and variance 
1 we deduce from Theorem 3 the following Corollary. 
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Corollary 1. Let H = H{N) be a function of N such that 

H XogH ^ 

■ +00, 



log N log TV 

as N — > 00. Assume that the hypothesis (20) holds for arbitrarily large K. Then the 

distribution of ip(n + H) — 'ip{n) — H for n < N is approximately normal with mean 
and variance HlogN, in the sense that the proportion of n < N for which 'i(j{n + 
H) — il){n) — H < cH log N tends to $(c) as N ^ 00, uniformly for \c\ < C. Here 
$(u) = /-oo ^"'^ cumulative distribution function of a normal random 

variable with mean and variance 1. 

For larger H, Theorem 3 furnishes only a limited number of moments and we cannot 
deduce a distribution result. However we expect that the contributions of the Ek{x;T>) 
cancel sufficiently so as not to overwhelm the main term: 

Conjecture 1. For each positive integer K , 



Mk{N; H) = {fiK + o{l))Ni^H\og - 



uniformly for {logNy+^ <H < N^-^ . 

This implies the weaker 

Conjecture 2. Suppose that (log A'")-'^"'"'' < H < N^~^ . The distribution ofip{x+H)—ip{x) 
for < X < N is approximately normal with mean H and variance HlogN/ H . 

Certainly Conjecture 2 does not hold when H ^ N, but perhaps it holds whenever 
H = o{N). It would be interesting to investigate more thoroughly what happens in this 
range. 

Hardy & Littlewood [10] provided heuristics that point toward the quantitative prime 
/c-tuple conjecture (1). In §4 we argue in the same spirit to obtain indications in favor of 
Conjecture 1. 

To obtain further support for our conjectures, we interpret the situation in terms of the 
zeros of the Riemann zeta function. We recall that Goldston & Montgomery [8] showed 
that if RH is true, then the stronger form {F{a) ~ 1) of the Pair Correlation Conjecture as 
formulated by Montgomery [13] is equivalent to the case K = 2 of the Conjecture above. 
In the same spirit, Chan [3] has shown (assuming RH) that Conjecture 1 is equivalent to 
the assertion that 



(22) / ( J2 cos(7loga;))''da;=(//fe + o(l))x(^logT)'^' 

0<7<T ^ 



Viewed in this way, we see that the Pair Correlation Conjecture asserts that the variance 
of the sum 

cos(7loga;) 

0<7<T 
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is the same as it would be if it were a sum of uncorrelated random variables, and Conjec- 
ture 1 asserts that this same sum has the same normal distribution that it would have if 
the terms were independent random variables. In somewhat the same vein, Bogomolny & 
Keating [1] used Hardy-Littlewood conjectures concerning primes to arrive at the n level 
correlation function of zeros of the zeta function. 

Freeman Dyson observed that the Pair Correlation Conjecture is analogous to known 
properties of the spacings of the eigenvalues of certain families of random matrices. We 
note that (22) has a similar analogue in random matrix theory. Let U{N) denote the 
classical compact group of unitary N x N matrices. For A G U{N), let e(6'i), . . . ,6(6^) 
denote the eigenvalues of A. Rains [19] has observed that if M is an integer, \M\ > N, 
then the point (M6'i, . . . , MOn) is exactly uniformly distributed in as A varies with 
respect to the Haar measure dA on U{N). It follows in particular that the distribution of 

N 

Re Trace = ^ cos(27rM^^) 

n=l 

is exactly the same as the distribution of 

N 

cos 2TTXn 

n=l 

where where the are independent random variables, each one uniformly distributed 
on [0,1]. It follows by the Central Limit Theorem that this distribution tends toward a 
normal distribution with mean and variance N/2. By an easy calculation it can also be 
shown that the k^^ moment of this distribution is if is odd, and is ~ //fc(A^/2)'^/^ if k 
is even. 

As for numerical studies. Brent [2] has compiled evidence not only for (1) but also for 
the stronger hypothesis (20). Odlyzko [18] and Forrester & Odlyzko [5] have found that 
the local distribution of the zeros of the zeta function fits well with predictions based on 
random matrix theory. The authors [15] have reported on numerical evidence in support 
of the conjectures. Finally, Chan [3; pp. 36, 49, 63] has assembled evidence in favor of (22). 

Cramer's model suggests that 

(23) n{x + (log x)") - n{x) ~ (log x)"-^ 

as a; — > oo with a fixed, a > 2. This, however, is known to be false, since Maier [11] showed 
that 

— 7r(a; + (logx)") -7r(a;) ^ ^ 
^-z^ (logo;)"-! < 

for any fixed a > (for general results of this nature see Granville & Soundararajan [7]). 
Presumably (23) is valid for most x, and the exceptions discovered by Maier are quite 
rare. Indeed Selberg [18] showed that on RH, (23) holds if a > 2 for almost all x, and 
Corollary 1 shows on hypothesis (20) that (23) holds if a > 1 for almost all x. As for 
longer intervals, suppose that a is fixed, < ck < 1. Cramer's model would predict that 



8 



MONTGOMERY & SOUNDARARAJAN 



'iIj{x + x"') — 'iIj{x) — x"^ is approximately normally distributed with mean and variance 
X°'\ogX as X runs over the range X < x < 2X. Our Conjecture 1 predicts normal 
distribution, but with a variance that is smaller by a factor of 1 — ct. Thus it seems that 
in this range, Cramer's model is not just occasionally inaccurate, but instead is actually 
inaccurate on average. 

1. Proof of Theorem 1 

Montgomery & Vaughan [16] devised a useful basic inequality (their Lemma 1), which we 
now quote. 

Lemma 1. Let ri,... ,rk be squarefree integers, set r = [ri, . . . ,rk], and suppose that 
any prime dividing r divides at least two of the r^. Then for any complex-valued functions 
Gi, . . . , Gfc defined on (0, 1] we have 

k 

E t[Gi{h/r,) 

bi,...,6fc i=l 

l<bi<ri 



Montgomery & Vaughan [17] have derived several variants of the above; an exposition 
of such variants is found in Chapter 8 of Montgomery [14]. For our present purposes a 
different type of variant is useful. 

Lemma 2. Let qi, ■ ■ ■ ,qk be squarefree integers, each one strictly greater than 1, and put 
d — [qi, ... ,qk\. Let G be a complex-valued function defined on (0, 1), and suppose that 
Gq is a nondecreasing function on the positive integers such that 



(24) 



Y,\Gia/q)\'<qGoiq) 



for all squarefree integers q > 1. Then 

k 



E IlG{ai/qi) <lflqiGo{qi)'/' 



0<ai<qi 



Proof. We write q^ — riSi where the Sj are pairwise relatively prime and any prime dividing 
[ri, . . . , rfe] divides at least two of the r^. That is, rj — {qi, Ilj^i Clearly d = rsi ■ ■ ■ s^ 
where r = [ri, . . . ,rfc]. The condition ^ai/qi e Z forces Sjjaj for all i. Hence, on writing 
ai = Sibi, we find that the left hand side above is 



E UGibi/n) 



bi,...,bk i- 

0<bi<ri 
E bi/nEZ 
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If there is an i for which = 1, then the conditions in the above sum cannot be fulfilled, 
the sum is empty, and there is nothing to prove. Thus we may assume that > 1 for all i. 
In Lemma 1 we take Gi{x) = G{x) for < a; < 1, and G{1) = 0. Thus by Lemma 1 and 
the hypothesis (24), the above is 

£ J n ('■■ E •■AG(b,/n)\f" < 1 n {r^G,{n)f\ 

i=l bi = l i=l 

Since G is nondecreasing, the above is 

k 



^^(^'^o(grf/^ = ^^(^.'^o(g.))'/^ 



< 

r 

i=l 1=1 



as desired. 



We now begin the main body of the proof of Theorem 1. We take k to be fixed, so 
that the dependence of implicit constants on k is suppressed. If k is odd, then the desired 
estimate is already found in (18) of Montgomery &; Vaughan [16]. Thus we may assume 
that k is even. From (12) and (14) it is clear that 

<li,---,qk ^1=1^^^*''/ ai,...,afc j=i 
Kgilg l<ai<qi 
{ai,qi) = l 
E ai/qi^'Z 

In Lemmas 7 and 8 of Montgomery & Vaughan [16], it is shown that all contributions to 
the above are 

\(l){q)J 

except for those terms for which the qi are equal in pairs, with no further equalities among 
the qi. There are (A; — 1)(A; — 3) • • • 3 • 1 = /Xfe ways in which this pairing can occur. Take 
the pairing to be qi = qk/2+ij ^ind set bi = ai + ak/2+i- Thus the terms that remain to be 
estimated are precisely 

fc/2 2 

(25) « E E Xl-'V'.i.) 

Kgik l<h,<q, 
qi distinct ^ bi/qiE'L 

where 

(26) JiP,r)= t E(^eC-^). 

o=l 
(a,r) = l 
(6-o,r) = l 
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First we show that the condition that the qi should be distinct in (25) can be dropped. To 
see this, put F{a) = min(/i, 1/||q;||) where ||^|| = min„gz \0 — n\ is the distance from 9 to 
the nearest integer. Thus 

(27) \E{a)\<F{a) 

for all a. Let Q denote the set of those /c-tuples {qi, . . . , q^) such that 1 < g^lg for all i, and 
with the property that among the qi there are three or more of them that are equal. In 
proving their Lemma 8 (see the treatment of T3), Montgomery & Vaughan [16] establish 
that 

qgQ \i=l '''Vyiy / ai,...,afe j=i V^V^/ 

l<ai<qi 
{ai,qi) = l 

Since this majorizes the difference between (25) and 

k/2 2 '^/^ 

E bi/qieZ 

it follows that we can continue with the above expression. Suppose that < bi < qi for 
exactly j values of i, and that bi = qi for the remaining k/2 — j values of i. Since there 
are ( '^^^ ) ways of choosing the j indicies, we see that the above is 

k/2 

(29) /^fc E ( )^2(g; h)''/^-m,{q; h) 



where Woi^q] h) = 1 and 



(30) wM.h)= E nS!% E \{Jip^^^^)- 

qi,---,qj i=l ^^y*^ bi,...,bj i=l 
l<9il9 0<6i<<ji 
E bi/qiEZ 

Here the term j = gives the desired main term. Thus it remains to show that the other 
terms are smaller. 

To prepare for an application of Lemma 2, we estimate J{b,r). By (27) we see that if 
< 6 < r/2 and r < h, then 

171 T + / IT TV + / ^ 

bib -a) ^ bia-b) ^ 

b/2<a<b ^ ' b<a<Zb/2 ^ ' 36/2<a<2r/3 

2 

(31) <^ log 26. 
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Here half the ranges of a have been omitted, since by symmetry they contribute the same 
amount as the hsted sums. Similarly, if < 6 < r/Zi and r > h, then 

(32) J{b,r)<^ I] + Yl —^rh. 

0<a<2r/h 2r/h<a<r/2 

Finally, ii r/h < b < r/2 and r > h, then 

^ ^ ^ hih-a) ^ a ^ a 

b/2<a<b-r/h ^ ' b-r/h<a<b b<a<b+r/h 

J.2 ^2 



b(a — b) ^ 0? 

b+r/h<a<4b/3 ^ ' 46/3<a<2r/3 



^2 

(33) <ylog(26VO- 
From (31) we see that if r < /i, then 

^ J(6,r)^«r^ 

0<6<r 

and from (32) and (33) we see that if r > /i, then 

0<b<r 

Altogether, 

(34) ^ J(6, r)^ -C r^min(r, /i). 

0<6<r 

On taking Go{r) — Chr"^ in Lemma 2, we find that 

0<6{<gi 

and hence 

d|g r|d p\q 

To apply this in (29), we need also a bound for V2{q'i h). To this end we note that 

d\q ^ ^ a=l d\q ' 

=''n(i+(^)«'"^- 
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(By a different method it can be shown that ^2(5; h) < hq/4>{q). See Hausman & Shapiro 
[11] and (3) of Montgomery & Vaughan [16].) From (30) we see that Wi{q;h) = 0, since 
the inner sum is empty. On applying the above estimates for 2 < j < k/2, we see that the 
expression (29) is 

f^,V,{q;h)'/' + 0{h^/'-\q/<l>{q)fy 



Here the error term is majorized by that in (15), so the proof is complete. 

2. Proof of Theorem 2 

We begin with two lemmas. 
Lemma 3. (Hardy-Littlewood) Let 

k 



ai,... ,afe i—i 
l<ai<qi 
(ai,qi) = l 



A{qi,... ,qk)= J2 KSZ^)- 

^ Hi ' 



If Qi = (M ^ith {UQi,UQi) = 1, then 

(35) A{q^, ...,qk) = A{q[, q'k)A{q'l, q^). 

For any prime number p, 

(-) E n^.(.,...,.)=(i-i)-'(i-^) 

qi,--- ,qk i=l ^V'**/ ^ ^ 

<li\p 

where ^'p(I') is the number of distinct residue classes modulo p found among the members 
of D = {di, . . . , dk}. Finally, 



(37) E n^i^(^i'---'^'^)i<-- 

l<qi<oo 



Proof. We follow the argument of Hardy & Littlewood [10, pp. 56-61], but with some 
helpful amplifications. We write 



a' a'/ 



^^^ + -L (modi). 

Qi Qi Qi 

By the Chinese Remainder Theorem, each reduced residue modulo qi corresponds to a 
pair a'-,a'^ of reduced residues modulo q'^ and q'/ , respectively. Also, ^Oj/gj £ Z if and 
only if e Z and J2(^i/Qi ^ ^- This gives (35). 
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If each Qi is either 1 or p, then 



^ ^ ^ gi >' I otherwise. 



r=l 1=1 

Thus the left hand side of (36) is 

■^r=igiv,gfc i=l ^^^'^ ai,...,afc j=i ^* 

l<ai<gi 
(ai,(7i) = l 

^ r=l i=l 0<o<p 

Here the innermost sum is p — 1 or —1, according as r = (mod p), or not. Thus if r = 
(mod p), then this factor of the product is 0. There are i^p{D) such values of r. For the 
remaining p — J>'p(2)) values of r, each factor of the product is p/ {p — 1). Hence the above 
is 



p-UpjV) / p \fc 

p — 1 / ' 



which gives (36). 

From (35) we see that 



Put D = Y[i<jidj ~ ^i)- -^o^ primes p\D we make no attempt to simplify the above 
expression. However, there are only finitely many such primes, and for p \ D, the di 
are distinct modulo p. For such primes we evaluate the factor more explicitly. Let d C 
{1, . . . ,k} with j = card J, suppose that qi = p for i & d, Qi = ^ for i ^ 3, and that the di 
are distinct modulo p. Then 

^fe.....*) = iE E i^)-l±m E K^))- 

^ r=l a.1,... ,ak r=l i=l l<ai<qi 

l<«i<9i (ai,gi)=l 

If i ^ d, then the innermost sum is 1 for all r. If z e d, then the innermost sum is p — 1 
if r = di (mod p), and —1 otherwise. Thus there are j values of r for which one factor is 
p — 1, j — 1 factors are —1, and all other factors are 1. For the remaining p — j values of 
r, there are j factors that are —1 and the remaining factors are 1. Thus the above is 

= -{j{p - i)i-iy-' + {p-j){-iy) = {-iy-\j - 1). 
p 
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Hence |A(gi, ... ,qk)\ = \ j — l\, so it follows that the expression (38) is 

n( E n^i-^te.-.*)i)xn(i+i:C)(^ 

p\Q ^qi,--- ,qk ^ p\Q \ j=2 ' 

Since this last product converges when extended over all primes, we have (37), and the 
proof is complete. 

From (35) and (36) we see that 

Qi,... ,qk i=l pig f 

qi\Q 



for any positive integer Q. By (37) it follows that 

k 



e(f)=i™ E n 



which by (39) is 

„„n(i-V-(i-^^). 



Thus the expressions (2) and (3) are equal. 

Suppose that 1 < < /i for all i. Then VpiT^) = k for all primes p > h, and thus if 
y > h, then 

(-) n(-r(-^)=n(i^o^(^))-+o*(j4)- 

P>y ^ ^ p>y ^ y toy 

Since I'piT)) > 1 for all p, it is evident that 

(41) n(i-V(i-^)«'(i°8'>)'-'- 

On combining this with (40), we see that 

(42) e(D) < (log/i)'=-^ 
From (5) it follows additionally that 

(43) 6o(a^)«(log/i)^-^ 
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From (39)-(41) we see that if 1 < di < h for all i and y > h, then 



By combining this with (5), we see also that 

(45) 6„(i,)= E ri^-<fe.-.*) + o4^^)- 

gi,--- ,qfe i=i y 

Lemma 4. Let E{a) be defined as in (14). If q is divisible by every prime number p < h? , 
then 

(46) E -<i'^d)h) = h'-hlogh + Bh + 0{h'/'+^) 

where B = 1 — Cq — log27r and Cq denotes Euler's constant. 
Proof. Since |£?(q;)P = '}2\m\<h(^ ~ |m|)e(mQ;), it follows that 



\E{a/d)\^^ {h-\m\)cd{m) 



a=l \m\<h 
(a,d) = l ' 



where Cd{m) is Ramanujan's sum. Now Cd(0) = 0(ci), and Cd{—m) = Cd{m), so the above 
is 

h 

^{d)h + 2 ^ (/i - m)cd{m). 
Hence the left hand side of (46) is 

m=l d\q ^ 

Here the sum over d is 

n(^+(^)=n(^%4T)n(-(^)- 

p\q p\q p\q 

p\m p\m 

Since q is divisible by every prime p < h^, the above is 

=n(i+^)n(i-(^)+o(v'>^)- 
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Here the main term is 



p\m p\m 

Goldston [7] has shown that 
h 

(47) 2 J] (/i - m)6({0, m}) = - hlogh + Bh + 0{h^'^+^), 

SO we have the stated estimate. 

It is worth noting that (47) can also be written in the form 

(48) &i{di,d2}) = h^ -h\ogh + Bh + 0{h^/^+^). 

di ,d2 
l<di<h 

The term d = 1 contributes h? — h to the left hand side of (46). Thus if q is divisible 
by every prime not exceeding /j,^, then 

(49) E -^^d)h) = -h\ogh + Ah^O{h'/'^^) 

d\q ' a=l 

d>l (a,d) = l 

where A = 2 — Cq — log 27r. 

We now begin the main body of the proof of Theorem 2. We apply (45) with y = h^^^ ^ 
and set Q = Y[p<y P- Thus 

(50) Rj,{h)= J2 t[jf^S{qu...,qk;h) + 0{1) 

qi,--- ,5fe i=i ^^^'>') 

Qi\Q 



where 



(51) S{qu...,qk;h)= Yl E KE 



ajdi 
Qi 

l<di'<h l<ai<5> 
di distinct {ai,qi)=l 
Yl a,i/qi€l, 



By comparing (50) with (12) we find that if the condition that the di should be distinct 
were omitted, then the main term in (50) would be exactly Vk{Q;h). The bulk of our 
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argument is devoted to an effort to remove this condition. Put Sij = 1 if cZj = dj, 5ij 
otherwise. Thus 

1 if the di are distinct; 



n (i-M=| 



otherwise. 



When the left hand side above is expanded, we obtain a hnear combination of products 
of the 5 symbols. Let A denote such a product, and |A| the number of factors in the 
product. We define an equivalence relation on these 5-products by setting Ai ~ A2 if Ai 
and A2 have the same value for all choices of di, . . . , (i^. For example, 5i,2^i,3 ~ ^1,2^2,3 ^ 
<^i,2<^i,3<^2,3- Given a partition (P = {Si, . . . , Sm} of the set {1, . . . , A;}, let 

M 



Ay = n n 



m=l i<j 



We see easily that every equivalence class of 5-products contains a unique Ay. Thus we 
have a bijcctivc correspondence between equivalence classes of 5-products and partitions 
of {1, . . . , /c}. For a partition (P, put 

(y)= E (-1)'^'- 



w ^ 

A~A 



Thus 



and it follows that 

M 

(52) 5(gi,...,5fc;/i) = EM3^) E H ^( E 



y ai,...,afe rn=l ie§ 

\<ai<qi 
{ai,qi) = l 
E ai/qiGl' 

where (P = {Si, . . . , Sm} 

If there is a prime p such that p\qi for exactly one i, then the condition ^ai/qi e Z 
cannot be fulfilled with (a^, g^) = 1, so the sum (52) is empty, and hence S{qi, . . . ,qk',h) = 
0. We therefore assume that each prime dividing [qi, . . . , qk] divides at least two of the qi. 

To facilitate our discussion of various types of partitions, we introduce some notation. 
Let M = {1, . . . , M}. For a partition = {Si, . . . , Sm}, put 

Ml = {m e M : cardSm = 1}, 
M2 = {m e M : cardS^ > 2}, 

^1 = U Sm, 

card S„i = l 



^2 = IJ Sm, 



mi 


= card Mi; 


777,2 


= cardM2; 


77l 


= card Ki ; 


772 


= card3Sf2. 



card S™ >2 
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Of course, mi + m2 = M, ni + n2 = k, and mi = ni. We first bound the contribution 
made by those partitions such that card Sm > 3 for some m. For m e Mi we use (27) to 
see that \E{ai/qi)\ < F{ai/qi). For m e M2 we use the trivial bound 



E 



ieSr. 



< h. 



Hence 



M 



ai,---,afe m=l ieS^ 

(ai,(ji) = l 



< h 



1712 



E U^i^ 



ai,...,afe ieXi 

l<ai<qi 
(ai,qi) = l 
E ailqi^I. 



For z G Ki we take Gi{x) = F{x) for < .x < 1, and G'i(l) = 0. For i e [N"2, we take 
Gi{x) — 1 for all X. Thus by Lemma 1 we see that the above is 



(53) 

\i q <h^ then 
If g > /i, then 



< 



m2 



[91,- •• ,9fc] , 



ie3Vi a=l 



n(.E-(^)7"n(..Er- 



1/2 



i6>r2 a=l 



q-1 



ek:) « E (!) <<«^. 



Q<a<q/2 



q-1 



0<a<q/h 



q/h<a<q/2 



q-1 



EKj) « E '^^^ E (!)<<.- 

a=l ' " .... 

Thus in any case, 

(54) ^F(-)'«5min(g,/.). 

a=l ^ 

Thus the expression (53) is 

^ qi---Qk ^ni/2+m2 

The members of ^2 are partitioned into m2 sets, each one containing at least two 
members, and at least one containing 3 or more numbers. Thus 

77-2= ^ cardSj„>l + 2 ^ l = l + 2m2, 

m€M.2 m€M.2 



and hence 



ni 77.1+77,2 — 1 k — 1 

Y + "^^- 2 = ^ 
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We also observe that 



^ 4>(qi) J \qu... ,qk] ~ ^ (b(q)J Z^il^ p-J 

,--^3fc ^'^^^^ ^ ^yfcj s^vy; d\Q p\d ^ 



Thus we have shown that 



M 



y gfc i=l ^''^^^ ai,...,afe ^=1 i£S^ ^» 

f55) cardS^<2 KgilQ l<ai<gi 

^ ' (ai,gi) = l 

Suppose that the partition J" consists of j doublcton sets and k — 2j singleton sets. Since 
no other 5-product is equivalent to Ay, and |Ay| = j, so w{y) — (— I)-'- The number of 
such partitions is 

Since the qi are interchangeable, we multiply by the above factor, and restrict our attention 
to one such partition: doubletons {z, i+j} for 1 < z < j and singletons {i} for < i < k. 
Thus the main term in (55) is 



0<j<k/2 ^'^^^ ^■'^\.^,<1. t\^(^^^ a,t',a, t\ + l * 

KqilQ l<ai<qi 
{ai,qi) = l 
J2 ai/qiEZ 

For 1 < i < j, let bi and be defined by the relations 

- = - + ^ (modi), l<bi<u, {bi,r,) = l, 
n qi qi+j 



and put 



(Ci,dij = l 
^ + ^^1 (mod 1) 
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Then the sum over the Qi in (56) is 

j , k 

r 

(f>{Qi) 



E E n«(^) E E n ^K^:) 

.-■■i^'j bi,...,b-j i=l * Q2j+i,---,qk 02^+1 ,a, 



(6i,ri) = l (ai,(/i)=l 

E ai/qiGZ 

We now separate those i for which = 1 from those with rj > 1. Let £ denote the number 
of i for which > 1. Since there are (;^) ways of choosing the £ values of the i from 
{1, . . . , j}, the above is 

j 



(58) J2{^^)Hiiy-'M{i) 
where 



^=0 



*nQ= E E n«(7) E E n fM^(^)- 

ri,— ,rebi,...,bei=l * q2j+i,— ,qk a2j+i,...,ak i=2j+l^^^^^ ^* 

KrilQ l<6i<ri KQilQ l<ai<qi 

(bi,ri) = l (04,90=1 

E fli/Qi+E bi/ri€Z 

We note that M(0) = V2fc-i(<5; /i)- 

Next we show that the contributions of £ > can be absorbed in the error term. If 
there is a prime p that divides exactly one of the numbers ri, . . . , r^, q'2j+i, • ■ ■ ^Qk, then the 
condition that 'Y^'^i/ Q.i '^'Yli^i/'^i ^ ^ cannot be fulfilled with (a^, qi) — {hi, ri) = 1, so the 
sum is empty, and the sum over the qi and vanishes. Thus we may restrict our attention 
to those choices of qi and for which every prime divisor of cZ = [q2j+i, ... ,qk,ri,... , r^] 
divides at least two of these numbers. Hence by Lemma 1, 

(-)w)«E E ifc.Ei<)irri(^EC)r- 



ri,... ,r(i q2j+i,... ,qk i=l 6=1 
KnlQ KqilQ 



In order to assess the above, we estimate H{h/r). We note that H{h/r) = if r f Q. 
Thus we suppose that r\Q, and that r > 1. For i = 1,2 we write di = Siti where Sjjr, 
U\Q/f- By the Chinese Remainder Theorem there exist unique Cj (mod Sj) and /j (mod 
ti) such that 

— + y = j- (mod 1), (ej, Sj) = (/j, U) = 1. 

Si li Ui 

From the conditions Ci/di + Cijd2 = hjr (mod 1), (6, r) = 1 it follows that 

^ = - (modi), [si,s2]=r, h ^ t2, fi + f2 = (mod ti). 

Sl S2 
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Put t = ti — t2- Hence 

Sl,S2 61,62 -l/'f'V i|Q/^ 

[si,S2]=r l<ei<Si 
(ei,Si) = l 
£1 +£2 = b (inod 1) 

For given si, S2, the number of pairs ei, 62 with the required properties is < </>((si, S2)) 
0(si)0(s2)/0(r). Hence 

[si,S2j=r 

<F - log/i. 

Vr / r 



By (54) it follows that 



r-l 



6^ 



2 /l. 



b=l 

On inserting this and (54) in (59), we find that 

E E in('.''^3"(-)iogft) n ^ 

ri\Q qi\Q 

« ""-^'"'"^ (">8 E K E 3"'-')' ( E 7^)'"' 

d\Q r\d q\d 
A k—2j-{-£ 

< h^'^-^^+'^/^ioghYH (1 + ) 

P\Q ^ 

< /j(fc-2i+^)/2+e_ 

From (57) we see that 

d>l 



Hence the expression (58) is 

(''E^)'^'-.(*'') + o('."-"'^+'). 



d>l 
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On inserting this in (56), we find that 

(60) «.w= E (2*)ff(-^E^)V.-..(Q;'» + (^<'-'"-'). 

0<i<fc/2 d\Q ' 

d>l 

We are at last prepared to make our appeal to Theorem 1. If A; is odd then k — 2j is 
odd, so there is no main term. Suppose that k is even. Then the main term is 

d>l 

d>l 

d>l 

by the binomial theorem. 

By (49) we see that the above is 

lik{-h\ogh + Ahf'^ + 0{h^^-^^/^+') . 

This gives the stated result. 

3. Proof of Theorem 3 

Clearly Mq[N\ H) — N . Thus we have the case K — Q oi (21) unconditionally, and with 
no error term. It is also convenient to dispose of the case K =1 before proceeding to the 
main argument. Since &{{h}) = 1, it follows from our hypothesis (20) that X]^=i -^('^ + 
h) = N + 0(7VV2+e). Hence En=i ^o(n + h) ^ TV^s+e^ thus by (19) we see that 
Mi{N;H) < HN^/^+^, which suffices. 

From now on we assume that K is fixed, K > 2, and we ignore possible dependence 
of implicit contstants on K. Let D — {di, . . . , dfc} be a set of /c distinct integers with 
1 < di < H for 1 < i < k. Suppose that the hi in (19) take the values di with multiplicities 
Mi. Then the right hand side of (19) is 

K JV' 1 N k 

(«) =E E (m.. mJm i: EnAo(n+*)- 

fc=l Mi,...,Mfe ^ di,...,dk "=1 

Mi>l l<d,<H 
E Mi=K di distinct 

Here the l/k\ is necessary because any permutation of (ii, . . . ,dk gives rise to the same set 
D. For positive integers m we put A^mip) — A(n)"^Ao(n). If M > 1, then by the binomial 
theorem 

M-l A/f-1 

Ao(n)^ = Ao(n)(A(n) - 1)^-^ = J] (_i)M-m-i /^^^ ')A^(n). 

m=0 
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On inserting this in (61), we find that 

k=l Mt,...,Mk 

Mi>l 

(62) E Mi=K 



mi,... ,mfe i=i •■ « 



0<mi<Mi 

where 

N k 

(63) Lfc(m)= J] J2 l[Am,{n + di). 

l<di<H 
di distinct 

To estimate the Lf^^m), we must distinguish between those i for which rrii — and those 
for which > 0. To this end we set 3C = {1, ... A;}, and introduce the following notation: 

!K = {i e % : rrii >1}, h = card IK; 

3 = {i e % : rrii = 0}, k - h = card J; 

dQOC, j = carda. 

Thus 

Yl Ao(n + c^i) n + ^i) = n ^o(^ + n ^^^in + di) + 1) 

3 i€3 
JC0C3C 

Prom our hypothesis (20) it follows that 

n<x ieJ ieM g 

ocgcx 

uniformly for < x < where = {di : i G 3}- With m = (mi, . . . , m^) fixed for the 
moment, write the above briefiy as f{x) = cx + 0(^N^/'^~^'^). Then 

N 



(64) E ( n + '^^)) ( n + diK^'^&in + di)r^-\log{n + di) - 1)) 

n=l ieO ie^K 

= 11 (log(x + di))^'-\\og{x + di) - 1) df{x), 



24 



MONTGOMERY & SOUNDARARAJAN 



which by integration by parts is 

/N 
Yl (log(x + di))"''-\\og{x + di) ~l)dx + 0(7VV2+s). 
- ie?{ 

For m > 0, A^(n) is nonzero only when n is a primepower, and 

Am{n) = A{n){lognr-\logn- 1) 

if n is prime. Thus if n is an integer such that n + di is prime for all i & !K, then the 
summand in (64) is 

k 

Those n for which n + di is a, higher power of a prime for one or more z e contribute an 
amount < ArV2+e to the sum (64). Thus the sum (64) is 

JV k 

(66) =J2 n^-^(^ + ^^) + 0{N'/^+'). 

n=l i=l 

Next we approximate the integral in (65) by a similar integral that is independent of the 
di- First we note that ii k = K, then Mj = 1 for all i, and hence = for all i, so that 
h = 0. Thus we may suppose that k < K. It is useful to note that 

k k 

(67) J] = J] m, < J](M, - 1) = K - k . 

ie'K i=l i=l 

If a; > 1 and 1 < d < iy, then 

log(a; + d) = logx + 0{d/x) = \ogx + 0{H/x). 
Thus the integrand is 

JJ ((loga;)"^'-\loga; - 1)) + 0{Hx-^{\ogN)^-^-^), 

and so the integral is Im{N) + 0(^H{logN)^~''^ where 

-AT 



On assembling our estimates, we find that 

N k 



Im{N) = / n (aogx)"^^-l(logX - 1)) dx. 



EnA^,(n + d,) = ( E &o{'^3))ilm{N)+0{H{\ogNf-'^))+0{N'/'+^). 
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By (43) we see that the first error term is < H{logN)^ < N'^/^+e ^ ArV2+e^ gince 
H < N^'^ and K >2. 

On summing both sides of the above over all choices of distinct di, we find that 

Lu{m)^Im{N) J2 E 6oPa)+0(if'=ivW^). 

3 di,... ,dk 
3C0C3C l<di<H 
di distinct 

Once the di have been chosen for i G ^J, there are {H — j){H — j — 1) ■ ■ ■ [H — k + 1) ways 
of choosing the remaining di. Hence the above is 

= Im{N) R,{H){H-j)...{H-k + l) + 0{H^N^'^+^). 

3 

The product in the sum is H''-^ + 0{H^-^-'^). By Theorem 2 we know that Rj{H) < 
{H log Hy/'^, and using this we see that (isolating the term J = J) 

J2 Rj{H){H-j)...{H-k + l) 

3 

3C3CX 

= Rk-h{H){H^ + 0{H^-^)) + 0( {H\ogHy/^H^-^) 

k-h+l<j<k 

= Rk-h{H)H^ + 0{{H log Hy^-^+^^/^H^-^). 

Therefore 

(68) Lkim) = Irn{N)(^Rk-h{H)Hf' + 0{iH\ogHy''-^+^y^H'^-^)') +0{H''N^/^+'). 

We insert the above in (62). In assessing the sizes of the various terms, it is useful to 
note that 

k 

(69) K = ^Mi = ^Mi + > 2card5{ + cardJ = 2/i+ (/c-Zi) = /i + Zc. 

i=i ieM iea 

By (67) we see that 

(70) ImiN) ~ N(\ogN)^ieo{^i < 7V(logiV)^-'=. 

First we show that terms for which h + k < K contribute a negligible amount to (62). 
Since Rk-h{H) < (if logif)^'^-'^)/^ by Theorem 2, it follows from (68), (69), and (70) that 

Lfc(m) < Ar(log7V)^-'=(iy logiy)^^-'^^/^^'^ ^ 
^ ^ ^ KlogNJ KlogNJ 



e 
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Thus the contribution of these terms to (62) is 

(71) « Ar(log AT)^ (_^)^''"'^^' + iy^ArV2+s_ 

Finally we consider those terms in (62) for which h + k = K. Since h < k, it follows 
that k > K/2. We also have card^K = h = K ~ k, and card J ~k — h = 2k — K. Since 
equality holds in (69), it follows that Mi = 2 for alH e 3-C and that Mj = 1 for all i e D. 
Thus rrii = Mi — 1 for all i, and for such m we have 

Im{N) = {logX - l)^-'=dx = lK-k{N), 

say. Hence 

Lfe(m) = lK-k{N)R2k-K{H)H''-^ + 0(iV(logiV)^-^(if logif)(2'=-^+i)/2ff^-'=-i) 

+ 0(iy^iV^/2+£) 

= lK-k{N)R2k-K{H)H^-^ + o(Ar(log7V)^(^)^''"'^^') + 0(i?^7VV2+s). 

Once /c is selected, there are precisely ( ) ways of choosing the set CK, and hence using 
(71) and the above, 

^72) K/2<k<K 

Suppose that K is odd. Then so also is 2k — K, and hence by Theorem 2 the main 
terms in (72) are 

«Ar(logAr)^/2iy^/2-V(7i^)+^ 

Thus in this case 

Mk{N; H) « iV(logiV)^/2iy^/2-V(7i^)+= + N{\ogN)^(^-^Y~^^'^ + H"" N^/^+' 

« Ar(fogAr)^/2iy^/2r_^r'/^'''^ + iy^ArV2+.. 

V log A* / 

Suppose that K is even. By Theorem 2 it follows that 



Mk{N;H) = h^/^ ^i^L^_J/^2.-x^^-.(iv)(-iogiy + ^)^-^/^ 

+ O (iV(log N)^'^H^/^ {\^T'^^''^ + if^iV^/^+^) . 



3g- 

Since /ifc = /c!/((/c/2)!2'^/^) when k is even, on writing k = K/2 + £ we see that the main 
term above is 

K/2 

P^kB^" E (^f)lKi2-m{- logi? + Ay. 
On taking the sum inside the integral, we obtain the stated result by the binomial theorem. 
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4. Heuristics in the manner of Hardy & Littlewood 

The reasoning here is conducted in the manner that Hardy & Littlewood [10] used to 
formulate their quantitative version of the prime A;-tuple hypothesis. Let 

N 

S{a) — A{n)e{na). 

n=l 

By the prime number theorem for arithmetic progressions we know that if (a, g) = 1, 
then S{a/q) ~ iJ,{q)N/(f){q) provided that q is not too large as a function of A^. By 
partial summation it follows that S{a) ~ ii{q)M{a — a/q)/(f){q) for a near a/q, where 

^iP) = '^n=i ^(''^f^) • -P^^ E{a) = Yl'^^i e{ma) , as in the proof of Theorem 1. If h is 
small compared with N, then S{a)E{—a) is approximately 

N h 

^ A(m + n)je{na) . 

n=l m=l 

Moreover, when a is near a/q, E{—a) is approximately E{—a/q). Thus when o: is a number 
for which the expression above is large, we expect that that it is approximately 

t. E{-a/q)M{a-a/q). 

l<q<N a=l 
{a,q)=l 

We subtract the contribution of the term a = q=l from both sides to see that 

N h M 

A(m + n)-/i)e(na)= ^ TTT E E{a/q)M{a - a/q). 

n=l m=l l<q<N a=l 

(a,q) = l 

Let F{a) denote the left hand side above. The A;-fold convolution of F with itself is 

N h , „ „ fc 

f A(m + n) — ^ j e{na.) — ' ' ' Y\. ^{(^i) dcni ■ ■ ■ da^ . 

}__, oti=ot (mod 1 j 

We set a = 0, and follow Hardy & Littlewood in assuming that the main term arises by 
the alignment of the peaks in the multiple integral on the right hand side. Thus we expect 
that 

N h k k 

n=l m=l qi,...,qk i=l^^^'^' ai,... ,ak i=\ 

l<qi<N l<ai<qi 
{ai,qi) = l 
E ai/<li€l 
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We note that \E{(3)\ x /i if \\/3\\ < 1/h, and that E{(3) < l/\\(3\\ if \\(3\\ > 1/h. The 
asymptotic size of the right hand side above could be determined by using the techniques 
used to prove Theorem 1 . At this point we are content to argue more informally. If k is odd 
then the terms do not make a very significant contribution. On the other hand, when k is 
even, we find 'diagonal terms' in which the qi are equal in pairs, with the corresponding 
tti being the negatives of each other. The pairings can be made in 



ways, so the contribution of these terms is 




If there are further equalities among the qi beyond this pairing, then the combinatorics must 
be adjusted, but such configurations contribute a lesser amount. Likewise, that the non- 
diagonal terms contribute a lesser amount can be demonstrated by using the techniques we 
used to prove Theorem 1. If q < h then the inner sum above is <^ Z^o<a<g ll*^/?!!"^ ^ 
but if q > h then the inner sum is approximately h(j){q). Since 



i:ff=ios.+o(i), 

<1<V 

it follows that the expression to be estimated is 

N\k/2 



likN[h\og ^ 



which supports the Conjecture. 
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